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INTRODUCTION 


Many  Navy  berths  that  are  located  in  estuaries  dominated  by  cohesive 
sediments  are  plagued  by  sedimentation.  Presently,  cutter-suction  dredges 
are  used  to  maintain  operable  depths  in  these  berths.  However,  the  in¬ 
convenience  of  these  dredges  and  the  ever-increasing  costs  associated 
with  dredging  and  dredge  spoil  disposal  have  prompted  the  Navy  to  con¬ 
sider  the  development  of  efficient  sediment  management  plans  for  each 
problem  harbor.  These  plans  would  include  sediment  resuspension  systems, 
sediment  exclusion  systems,  and  an  optimized  dredging  schedule.  To 
initiate  these  solutions,  it  is  necessary  to  have  a  rational  means  of 
analyzing  the  sediment  transport  problems  thac  currently  exist  in  each 
harbor  as  well  as  a  means  of  predicting  how  the  sediment  transport  might 
vary  with  future  harbor  developments. 

Most  available  predictive  models  of  estuarial  fine-sediment  trans¬ 
port  are  based  upon  theoretical  formulations  supported  by  laboratory 
flume  experiments.  Since  field  data  pertaining  to  these  transport  pro¬ 
cesses  are  particularly  scarce,  it  is  unknown  how  well  these  laboratory 
estimates  relate  to  the  actual  process. 

This  work  was  funded  through  the  Independent  Research  Program  at 
the  Naval  Civil  Engineering  Laboratory  (NCEL)  by  the  Office  of  Naval 
Research  (ONR) . 


OBJECTIVE 

The  commonly  used  models  for  erosion,  deposition,  and  transportation 
of  suspended  sediments  will  be  reviewed.  Many  of  these  models  require 
parameters  that  are  either  poorly  understood  or  difficult  to  quantify. 

The  objective  of  tl  is  study  is  to  investigate  the  true  character  of  these 
parameters  for  one  estuarine  environment,  and  to  use  this  information  as 
a  basis  from  which  to  comment  on  the  appropriateness  of  the  estimation 
methods  currently  available  for  cohesive  sediment  transport  modeling. 


APPROACH 

A  field  experiment  was  undertaken  to  gather  the  data  necessary  to 
address  the  quantitative  validity  of  the  most  popular  methods  for  pre¬ 
dicting  cohesive  sediment  transport. 


BACKGROUND 

There  are  three  processes  involved  in  estuarial  cohesive  sediment 
transport:  (1)  the  erosion  of  the  sediment  bed,  (2)  the  deposition  of 

suspended  sediment  onto  the  bed,  and  (3)  the  transportation  of  suspended 
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sediments  laterally  through  the  water  column.  Each  of  these  processes 
has  to  be  modeled  in  a  reasonable  fashion  in  order  to  predict  the 
shoaling  characteristics  of  a  harbor. 

Erosion 


The  process  of  erosion  involves  the  resuspension  of  material  from 
the  sediment  bed.  This  can  occur  either  as  surface  erosion,  where  one 
particle  at  a  time  is  removed,  or  as  bulk  erosion,  where  large  masses  of 
sediment  are  resuspended  at  once.  In  quasi-steady  flow  environments, 
surface  erosion  typically  dominates.  There  is  a  wide  diversity  of  equa¬ 
tions  that  have  been  used  to  describe  the  rate  of  surface  erosion  (Mehta, 
et  al.,  1982).  Of  these,  the  formulation  presented  by  Partheniades 
(1962)  is  the  most  common: 


E  = 


M 


T 

C 


where:  E 

M 


T 

C 


rate  of  surface  erosion 
erosion  rate  coefficient 
bed  shear  stress 
critical  shear  stress 


(1) 


Two  of  the  three  parameters  in  Equation  1,  the  erosion  rate  coeffi¬ 
cient  and  the  critical  shear  stress,  are  not  well  understood. 

A  few  attempts  have  been  made  to  quantify  the  erosion  rate  coeffi¬ 
cient,  M.  The  estimates  of  M  resulting  from  these  investigations  vary 
by  two  orders  of  magnitude.  Ariathurai  and  Arulanandan  (1978)  placed 
remolded  cylindrical  samples  of  differing  sediment  types  and  pore  fluid 
compositions  into  a  bath  of  eroding  fluid.  They  found  that  the  erosion 
rate  coefficient  varied  from  3  to  30  milligrams/square  centimeter-minute. 
The  value  of  the  coefficient  was  most  significantly  affected  by  the  cation 
exchange  capacity  (a  measure  of  the  soil  chemistry)  and  the  sodium  absorp¬ 
tion  ratio  (a  measure  of  the  water  chemistry).  More  recently  Ockenden 
and  Delo  (1988)  conducted  a  series  of  flume  experiments  to  determine  the 
erosion  rate  coefficient.  Using  an  annular  flume,  they  imposed  a  uniform 
bed  shear  stress  upon  a  naturally  deposited  sediment  bed.  Based  on  these 
experiments  they  suggest  that  the  erosion  rate  coefficient  is  a  function 
of  the  proportion  of  sand  within  the  sediment  (a  sieved  natural  sediment 
was  used).  Values  ranged  from  0.36  milligrams/square  centimeter-minute 
for  0  percent  sand  to  0.06  milligrams/square  centimeter-minute  for  11 
percent  sand. 

The  critical  shear  stress  is  a  rather  nebulous  term  that  addresses 
the  strength  properties  of  the  sediment  bed.  Parchure  and  Mehta  (1985) 
have  shown  that  this  value  is  not  constant,  but  varies  with  depth  into 
the  sediment  bed.  Based  on  laboratory  tests,  they  suggest  that  this 
variation  lends  itself  to  a  three-zoned  description  (see  Figure  1). 
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Each  of  these  zones  is  defined  by  a  reference  shear  strength:  the  upper 
zone  is  defined  by  the  shear  strength  at  the  sediment  bed  surface  (z  =  0), 
t  ,  which  is  the  critical  shear  stress  of  erosion  initiation;  the  middle 
zone  is  defined  by  the  shear  strength  at  the  point  of  greatest  curvature 
of  the  shear  strength  versus  depth  curve,  t  ;  and  the  lowest  zone  is 
marked  by  the  maximum  shear  strength,  t  .  S^he  range  of  values  obtained 
for  redeposited  clay  in  the  laboratory  for  each  of  these  shear  strengths 
is  listed  in  Table  1. 

Deposition 

The  process  of  deposition  involves  the  settling  of  particle  aggre¬ 
gates  from  the  water  column.  These  aggregates  adhere,  in  whole  or  in 
part,  to  the  sediment  bed.  Krone  (1962)  describes  this  process  as 
follows : 


dc 

dt 


■  -( 


cW  P 
s 


(2) 


where:  dc/dt  =  deposition  rate  (time  rate  of  change  of  the 

suspended  sediment  concentration) 

c  =  suspended  sediment  concentration 

Wg  =  aggregate  settling  velocity 

P  =  probability  of  a  particle  sticking  to  the  bed 

y  =  water  depth 

The  probability  of  a  particle  sticking  to  the  bed  is  defined  as: 


P  =  1 


- 

Tcd 


(3) 


where:  x  =  bed  shear  stress 

w 


t  ,  =  critical  shear  stress  of  deposition  (x  .  <  x  ) 
cd  r  cd  w 


As  with  the  erosion  process,  the  critical  shear  stress  parameter 
for  the  deposition  process  is  not  well  understood.  The  critical  shear 
stress  of  deposition  is  defined  as  the  shear  stress  below  which  all 
suspended  sediment  eventually  deposits.  Its  value  has  been  estimated  by 
conducting  laboratory  flume  experiments  using  sediments  sampled  from 
various  locations.  The  differences  i"  these  values  may  be  related  to 
grain  characteristics,  but  no  generally  accepted  relationship  is 
available.  Table  2  lists  these  values. 
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Transportation 


The  process  of  transportation  involves  the  lateral  movement  of  the 
suspended  sediment  load  through  the  water  column;  this  process  can 
alternately  be  defined  as  the  mass  flux  of  sediment  through  discrete 
control  volumes  within  the  estuary  (see  Figure  2).  The  mass  flux,  F, 
through  a  discrete  surface  is  defined  as: 

F  =  l  c  V  •  n  dS  (4) 


where:  c  =  suspended  sediment  concentration 


V  =  horizontal  fluid  velocity  vector 


n  =  unit  vector  normal  to  dS 


dS  =  incremental  area  on  the  surface,  S 


By  assuming  two-dimensional  flow  and  using  vertical  distributions 
for  both  the  suspended  sediment  concentration  and  the  horizontal  fluid 
velocity,  one  reduces  the  above  equation  to: 


■/ 


CV  dz 


where:  C  =  vertical  distribution  of  the  suspended  sediment 

concentration 

V  =  vertical  distribution  of  the  horizontal  fluid  velocity 
dz  =  incremental  depth 
h  =  total  depth  of  water 

The  vertical  distribution  of  the  suspended  sediment  concentration  is 
normally  defined  by  the  Rouse  equation  (Rouse,  1938)  as: 


(*?)  (sh) 


W  /kU. 
s  * 


where:  C  =  suspended  sediment  concentration  at  a  height,  z,  above 

the  sediment  bed 

C  =  reference  concentration  at  a  height,  a,  above  the 
sediment  bed 


h  =  total  depth 
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z  =  height  above  the  sediment  bed 

a  =  height  above  the  bed  of  the  reference  concentration 
=  settling  velocity 
k  =  Von  Kantian' s  constant 
U.  =  friction  velocity 

This  formulation  requires  that  a  reference  concentration  be 
specified  for  each  point  at  which  the  mass  flux  is  required.  Since  it 
would  not  be  practical  to  measure  these  values  for  a  predictive  model, 
some  attempts  at  developing  a  theoretical  basis  have  been  published 
(Adams  and  Weatherly,  1981;  Shi,  et  al.,  1985).  In  a  recent  field  study 
however,  Sternberg,  et  al.  (1986)  indicated  that  these  predictions  oi 
the  reference  concentration  were  not  accurate  enough  to  use  in  mass  flux 
estimates . 

The  other  parameter  in  Equation  6  that  is  not  well  understood  is 
the  settling  velocity.  In  most  modeling  it  is  taken  to  be  a  constant, 

much  as  three  orders  of  magnitude 
1988)  (see  Figure  3). 
fluid  velocity  is  defined  by  the 
flow  over  a  flat  plate: 

(7) 

z,  above  the  sediment  bed 


z  =  height  above  the  sediment  bed 
z^  =  bed  roughness  height 
The  friction  velocity  is  defined  as: 


yet  tests  show  that  it  can  vary  by  as 
within  the  same  estuary  (van  Leussen, 
The  vertical  distribution  of  the 
turbulent  boundary  layer  equation  for 


where:  U  =  mean  velocity  at  height, 

UA  ~  friction  velocity 
k  =  Von  Karman's  constant 


U 


* 


>r 


(8) 


where:  T  =  bed  shear  stress 

o 

p  =  fluid  dens i ty 


Von  Karman's  constant  for  a  homogeneous  fluid  is  0.4.  It  has  been 
shown  that  this  value  varies  in  flows  in  which  there  is  a  density  gradient 
(Vanoni  and  Nomicos,  1959).  It  is  not  certain  what  effect  a  suspended 
sediment  induced  density  gradient  would  have  on  this  constant.  Adams 


and  Weatherly  (1981)  and  Sternberg,  et  al.  (1986)  have  used  a  Richardson 
damping  scheme  to  simulate  the  effects  of  stratified  flow,  but  there  is 
little  experimental  verification  for  this  approach. 

FIELD  EXPERIMENT 

A  field  experiment  was  conducted  by  the  Naval  Civil  Engineering 
Laboratory  (NCEL)  from  23  to  25  September,  1987.  Suspended  sediment 
concentration  and  current  velocity  data  were  collected  in  the  bottom 
boundary  layer  of  an  estuarine  tidal  channel  over  a  duration  of  five 
tidal  cycles.  These  data  were  used  in  an  attempt  to  evaluate  the 
parameters  of  interest  in  the  above-mentioned  equations. 

Si  te 


The  field  study  was  conducted  in  Mare  Island  Strait,  an  estuarine 
tidal  channel  Lhat  acts  as  the  principal  effluence  of  the  Napa  River 
into  the  northeastern  corner  of  the  San  Francisco  Bay  complex.  The 
normal  tidal  flux  within  this  narrow  estuary  greatly  exceeds  the  fresh¬ 
water  effluence  throughout  most  of  the  year.  Tide  data  are  shown  in 
Figure  4.  The  surface  wave  activity  within  the  strait  is  negligible, 
and  the  tide-generated  current  is  the  predominant  driving  force  for 
sediment  motion.  The  depth  in  a  250-meler-wide  channel  and  in  numerous 
shipyard  berthing  areas  is  maintained  at  10.8  meters  below  the  mean 
lower  low-water  level  by  a  year-round  dredging  effort.  The  sediment  at 
Mare  Island  consists  of  46  percent  silt,  42  percent  clay,  and  12  percent 
sand  (Malloy,  1980).  Roughly  2.5  million  cubic  meters  of  this  sediment 
is  removed  annually. 

System 

A  sliding  instrumented  sting,  mounted  upon  a  bottom- landing  frame¬ 
work,  was  lowered  onto  a  flat  portion  of  the  artificially-deepened  channel 
approximately  70  feet  off  the  quay  wall  that  forms  the  western  boundary 
of  the  waterway  (Figure  5).  This  sting  was  oriented  perpendicular  to 
the  channel  axis  so  as  to  hold  the  instruments  in  an  environment  of  un¬ 
obstructed  flow  (Figure  6).  By  moving  the  sliding  sting  in  the  vertical 
plane,  via  a  hydraulic  cylinder,  it  was  possible  to  use  a  relatively 
small  number  of  instruments  to  profile  the  bottom  boundary  layer  of  the 
estuary.  Instruments  were  mounted  on  the  sting  and  used  as  follows: 

•  Two  Marsh  McBirney  electromagnetic  current  meters  (1-1/2-inch 
ball)  were  used  to  measure  the  horizontal  velocity  vectors  in 
the  water  column. 

•  Five  optical  backscatterance  sensors  were  used  to  determine  the 
suspended  sediment  concentration  in  the  water  column. 

•  A  Benthos  plankton  camera  was  used  to  determine  the  vertical 
distribution  of  floe  sizes  in  the  water  column.  Problems  were 
encountered  with  the  operation  of  the  device,  therefore  no  data 
were  obtained. 
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•  A  bottom  sensor,  based  on  the  design  presented  by  Krone  (1972), 
was  employed  to  reference  the  above-mentioned  measurements  to 
the  location  of  the  bottom. 

•  A  differential  pressure  transducer  was  used  to  measure  the 
excursion  of  the  sliding  sting  relative  to  its  highest  point  of 
excursion  (a  fixed  point  on  the  framework). 

Technique 

Vertical  profiles  of  the  fluid  and  sediment  parameters  were  taken 
at  half-hour  intervals  during  the  time  that  che  system  was  monitored 
manually.  The  profiling  was  performed  in  the  following  manner:  The 
sting  was  moved  downward  until  the  bottom  sensor  indicated  that  the  sedi¬ 
ment  bed  had  been  reached.  In  this  position  the  differential  pressure 
sensor  was  zeroed  and  the  first  set  of  readings  was  taken.  The  sting 
was  moved  up  to  the  next  station  and  a  new  set  of  readings  was  taken. 
Figure  7  shows  the  positions  of  the  stations  at  which  the  readings  were 
taken. 

Each  profile  (consisting  of  ten  stations)  required  about  3  minutes 
to  complete.  It  was  assumed  that  steady  fluid  flow  conditions  existed 
during  this  time.  To  insure  that  this  was  the  case,  readings  were  taken 
as  the  sting  moved  both  upward  and  downward.  The  results  were  compared, 
and  no  appreciable  offset  was  seen  in  any  of  the  profile  data  thereby 
suggesting  that  the  assumption  was  valid. 

In  addition  to  the  profile  data,  a  time  series  of  data  was  taken 
every  10  minutes  at  one  vertical  station  (with  the  bottom  sensor  located 
5  centimeters  above  the  bed)  for  the  duration  of  the  experiment. 


RESULTS 


Velocity  Data 

Current  velocity  data  were  collected  by  two,  two-axis  electromag¬ 
netic  current  meters  mounted  in  a  vertical  stack.  These  sensors  were 
calibrated  for  unidirectional  flow  using  a  tow  channel.  In  an  attempt 
to  reduce  any  data  scatter  that  might  have  been  caused  by  small-scale 
turbulence,  each  data  point  that  was  recorded  was  an  average  of  ten 
measured  points  taken  over  an  interval  of  3  seconds.  The  horizontal 
velocity  vector  for  each  position  was  calculated  as  follows: 


=  (u 


2,1/2 

v  ) 


0  =  tan 


(9) 


where:  |u| 
u 

v 

0 


magnitude  of  the  velocity  vector 

current  velocity  along  the  channel  axis 

current  velocity  perpendicular  to  the  channel  axis 

direction  of  the  velocity  vector  relative  to  the 
channel  axis 
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It  was  found  that  there  was  substantial  scatter  in  the  velocity 
data,  especially  near  slack  tides  where  the  magnitude  of  the  velocity 
vector  was  near  that  of  the  instrument  resoli  tion  (commonly  accepted  as 
roughly  3  centimeters/second).  This  scatter  was  probably  due  to  both 
the  insensitivity  of  the  instrument  relative  to  such  low  velocities  and 
to  fluctuations  caused  by  large-scale  turbulent  eddies.  Data  for  the 
measured  velocities  appear  in  Appendix  A. 

Thp  friction  velocity  for  each  data  set  was  comouted  by  using  the 
Von  Karman-Prandt 1  logarithmic  expression  given  in  Equation  7.  A  least- 
squares  linear  regression  fit  was  made  to  the  semilog  plot  of  velocity 
versus  the  log  of  the  height  above  the  bed.  An  example  profile  is  shown 
in  Figure.  8.  The  slope  of  this  line  is  k^/U,  and  the  intercept  is  z  . 

A  Von  Karman's  constant  of  0.4  was  used  in  the  analysis.  Due  to  the 
scatter  in  the  data,  it  was  deemed  inappropriate  to  attempt  to  look  at 
trends  beyond  the  linear.  This  rules  out  any  effort  to  describe  vari¬ 
ations  in  the  value  of  the  constant  due  to  density  stratification  (i.e., 
Richardson  damping). 

Friction  velocity  values  were  also  obtained  from  the  time  series 
data  in  a  similar  manner.  These  values,  obtained  from  only  two  data 
points  each,  were  compared  to  that  of  the  profiles  (20  to  25  points 
each)  to  determine  the  error  introduced  by  reducing  the  number  of  data 
points.  In  general,  it  was  found  that  the  two-point  estimates  were 
reasonably  reliable.  Data  for  all  calculated  friction  velocities  appear 
in  Appendix  B.  Based  upon  these  friction  velocities,  values  for  the  bed 
shear  stress  were  computed  using  Equation  R.  These  values  will  be  used 
in  conjunction  with  the  suspended  sediment  concentration  data  and  Equa¬ 
tion  1  to  determine  field  values  of  the  critical  shear  stress  and  erosion 
rate  coefficient. 

Concentration  Data 

Five  optical  backscatterance  sensors  were  used  to  measure  the 
suspended  sediment  concentration.  This  sensor  determines  the  amount  of 
backward  scattering  of  infrared  light  caused  by  a  water  sample  (Downing, 
et  al.,  1981),  and  relates  it  to  the  particulate  concentration  via  in¬ 
strument  calibration.  The  sensors  described  herein  were  calibrated  in 
the  laboratory  both  before  and  after  the  experiment  using  samples  of 
sediment  taken  from  the  sediment/water  interface  of  the  experiment  site. 
Figure  9  shows  an  example  calibration  for  these  sensors. 

Each  recorded  data  point  was  an  average  of  ten  measured  points 
taken  over  3  seconds.  This  approach  was  adopted  in  an  attempt  to  reduce 
any  scatter  that  might  have  been  caused  by  small-scale  turbulence.  Both 
profile  and  Lime  series  data  were,  collected.  These,  data  are  preserved 
in  Appendix  C. 

From  a  comparison  of  the  shear  stress  time  series  with  that  of  the 
suspended  sediment  concentration,  it  was  noted  that  reasonable  correla¬ 
tions  were  only  observed  for  the  ebb  tide  data.  The  flood  tide  concen¬ 
tration  events,  which  included  much  larger  suspended  sediment  concentra¬ 
tion  levels  than  did  the  ebb,  did  not  correspond  at  all  with  the  higher 
shear  stresses  as  should  be  expected.  Rased  upon  past  studies  of  the 
siltation  mechanism  at  Mare  Island,  the  following  explanation  is  offered. 
Mare  Island  Strait  is  a  small,  deep  channel  that  is  appended  to  Carquinez 
Strait  very  near  to  its  effluence  into  San  Pablo  Bay  (see  Figure  5). 


8 


Carquinez  Strait  is  a  deep  channel  with  extremely  fast  currents  (in  some 
areas  it  has  scoured  to  a  depth  of  70  feet).  San  Pablo  Bay  is  a  shallow 
muddy  embayment  of  substantial  fetch.  The  waves  generated  by  the  typical 
onshore  winds  suspend  mud  in  the  flats  near  the  mouth  of  the  Strait. 

The  combination  of  the  intense  mixing  in  Carquinez  Strait  and  the  wave¬ 
generated  mud  suspension  in  San  Pablo  Bay  creates  an  ample  source  of 
suspended  sediment  that  is  washed  into  Mare  Island  Strait  during  the 
flood  tide  (Krone,  1987).  Using  this  information,  it  is  argued  that  the 
flood  tide  concentration  events  were  not  resuspension  events  caused  by 
local  surface  erosion,  but  rather  an  advection  of  suspended  sediment 
that  was  resuspended  outside  the  zone  of  the  assumed  uniform  flow.  With 
this  argument  as  justification,  the  flood  tide  data  were  removed  from 
consideration  in  the  analysis  that  is  described  below. 

Because  suspended  sediment  concentrations  in  the  lower  regions  of 
the  bottom  boundary  layer  were  smaller  than  expected,  the  lower  two 
optical  backscatterance  sensors  were  out  of  their  operating  range  for 
much  of  the  time.  They  only  operated  during  resuspension  events.  Yet, 
it  is  necessary  to  have  these  values  to  obtain  the  total  suspended  load 
for  a  given  time.  Due  to  this  problem  it  was  necessary  to  attempt  to 
extrapolate  the  concentration  curves  to  cover  the  lower  region  of  the 
bottom  boundary  layer  for  times  of  low  suspended  sediment  concentrations. 
This  extrapolation  was  based  on  Equation  6.  A  least-squares  linear  re¬ 
gression  fit  was  made  to  the  log- log  plot  of  concentration  versus  dimen¬ 
sionless  depth,  [ (h-z)/z] [ a/(h-a) ] .  The  slope  of  this  line  is  equal  to 
the  exponent,  W  /kUVf.  Based  on  this  exponent,  the  concentrations  were 
found  for  the.  two  stations  in  which  data  were  missing. 

The  rate,  of  surface  erosion  has  dimensions  of  mass  eroded  per  unit 
time.  This  can  also  be  expressed  as  the  time  derivative  of  the  total 
suspended  sediment  load.  Thus,  by  substituting  this  expression  into 
Equation  1  and  rearranging: 
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(_V] 

d  (TSL) 
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dt 
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(10) 


,  d(TSL) 
where 

dt 


=  the  time  dei  ative  of  the  total  suspended  load. 


A  least-squares  linear  regression  fit  was  made  to  the  plot  of  bed 
shear  stress  versus  the  time  derivative  of  the  total  suspended  load. 

The  slope  of  the  line  is  T  /M,  and  the  y-intercept  is  the  critical  shear 
stress  (Figure  10).  Based  on  this  information,  the  critical  shear  stress 
(x  )  of  the  surface  layer  of  sediment  in  Mare  Island  Strait  is  1.44 
dynes/square^cent imeter .  The  erosion  rate  coefficient  (M)  is  estimated 
as  3.5  x  10  mi  1 1 igrams/square  cent imeter-minute. 

As  a  check,  the  critical  shear  stress  (x  )  was  estimated  in  a 
second  manner,  utilizing  data  for  a  concent raE i oi  increase  for  each  of 
the  five  measured  erosion  events.  The  bed  shear  stress  value  corres¬ 
ponding  to  the  inception  of  an  increase  in  suspended  sediment  concen¬ 
tration  was  noted.  Because  both  the  concentration  sensor  position  and 
the  finite  sampling  interval  allow  a  time  lag,  it  was  assumed  that  the 
increase  in  concentration  actually  began  a  short  time  before  the  sensors 
detected  it  (sometime  between  the  detect  ion  reading  and  the  previous 
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one).  Therefore,  the  bed  shear  stress  at  both  the  point  of  detection  of 
the  increase  in  concentration  and  at  the  previous  point  were  considered. 

The  critical  shear  stress,  as  defined  by  Equation  1,  is  a  constant  for  a 
given  sediment  type  and  depth  into  the  seabed.  Therefore,  it  is  reasonable 
to  assume  that  the  same  value  of  bed  shear  stress  that  is  critical  for 
this  given  situation  must  be  present  in  all  five  of  the  measured  erosion 
events.  Allowing  for  some  measurement  uncertainty,  the  only  range  of 
bed  shear  stresses  that  is  present  in  all  five  measured  erosion  events 
is  1.0  to  1.3  dynes/square  centimeter  (Figure  11).  This  is  in  reasonable 
agreement  with  the  more  precise  value  of  1.44,  found  with  the  previous 
method. 

The  other  key  parameter  in  Equation  10,  the  erosion  rate  coefficient 
(M),  is  probably  slightly  underestimated.  The  flux  of  sediment  into  the 
upper  water  column  was  not  taken  into  consideration  in  the  calculation 
of  the  total  suspended  sediment  load.  This  flux  would  occur  as  a  gradual 
process  through  the  later  stages  of  the  resuspension  event.  This  suggests 
that  the  time  derivative  of  the  total  suspended  load,  especially  in  the 
higher  portions  of  the  curve,  is  underestimated.  An  attempt  has  been 
made  to  estimate  the  magnitude  of  this  error  by  using  the  Rouse  equation 
to  extrapolate  the  sediment  concentration  data  into  this  portion  of  tne 
water  column.  These  efforts  suggest  that  the  error  is  up  to^lS  percent, 
thu^  offering  a  possible  range  of  values  of  M  from  3.5  x  10  to  4.0  x 
10  milligrams/square  centimeter-minute. 

Apparent  Settling  Velocity 

An  apparent  settling  velocity  was  estimated  using  Equation  6.  A 
least-squares  linear  regression  fit  was  made  to  the  log-log  plot  of 
concentration  versus  dimensionless  depth,  [ ( (h-z)/z; (a/(h-a) ) ] .  The 
slope  of  this  line  is  equal  to  the  exponent,  W  /kU^..  Both  k  and  are 
known  for  each  set  of  measurements;  therefore ,' the  settling  velocity,  W 
can  be  obtained.  This  procedure  was  performed  for  both  the  profile  data 
and  the  time  series  data.  Fair  agreement  was  found  between  tnem.  A 
graph  of  these  values  can  be  found  in  Appendix  D. 

The  apparent  settling  velocity  data  (profile  data  only)  were 
compared  to  concentration,  velocity,  and  friction  velocity  data  in  an 
attempt  to  correlate  the  apparent  settling  velocity  to  more  easily 
measured  parameters.  No  correlation  was  found  with  either  the  fluid 
velocity  or  the  friction  velocity.  A  weak  correlation  was  found  between 
the  apparent  settling  velocity  and  the  suspended  sediment  concentration. 
This  correlation  is  shown  in  Figure  12.  A  number  of  settling  velocity 
experiments  have  been  conducted  in  quiescent  waters  (Krone,  1962;  Owen, 
1970;  Kranck,  1986)  that  have  considered  similar  correlations.  Results 
show  that  the  sediment  concentration  is  a  strong  factor  in  determining 
the  settling  velocity,  but  that  the  relationship  is  not  linear.  The 
exponents  cited  have  ranged  from  1.0  to  2.5  f->r  saline  water  (van  Leussen, 
1988).  The  field  data  presented  herein  are  not  of  sufficiently  high 
resolution  to  determine  whether  or  not  a  better  fit  would  be  obtained  by 
using  an  exponent  other  than  one. 

In  standard  settling  velocity  experiments  in  quiescent  waters 
the  dominant  cause  of  flocculation  is  collision  of  particles  due  to  the 
different  settling  rates  of  different  sized  aggregates.  It  has  been 
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suggested  (Krone,  1972)  that  this  cause  is  actually  of  secondary  impor¬ 
tance  when  compared  to  the  collisions  caused  by  internal  shearing  of  the 
fluid  flow.  Therefore,  the  effects  of  the  friction  velocity  were  included 
in  the  apparent  settling  velocity  versus  sediment  concentration  comparison. 
A  much  stronger  correlation  was  found  than  was  seen  when  comparing  the 
apparent  settling  velocity  with  concentration  alone  (Figure  13).  This 
result  supports  Krone's  inference,  and  casts  a  shadow  of  doubt  over  all 
of  the  efforts  to  extrapolate  a  field  settling  velocity  from  experiments 
conducted  under  quiescent  conditions. 

Some  of  the  time  series  data  sets  had  five  points  available  from 
which  to  obtain  the  apparent  settling  velocity.  Most  had  only  three 
points.  An  apparent  settling  velocity  was  calculated  for  the  five-point 
data  sets  using  both  five  points  and  the  three  points  that  were  available 
in  the  rest  of  the  data  sets.  It  was  noted  that  the  apparent  settling 
velocity  values  were  generally  higher  when  five  points  were  used  in  the 
estimation  process.  This  suggests  that  the  settling  velocity  may  be 
larger  in  the  lower  portions  of  the  water  column.  This  inference  is 
supported  by  data  presented  in  Sternberg  et  al.  (1986)  that  indicate 
that  aggregate  particle  size  is  typically  "  irger  in  the  lower  portions 
of  the  water  column.  If  the  settling  ve.  ity  varies  significantly  with 
depth,  the  equations  describing  deposition  and  transportation  may  require 
revis ion . 


SUMMARY 

A  field  experiment  was  conducted  to  gather  data,  which  were  then 
used  to  evaluate  methods  for  predicting  cohesive  sediment  transport. 

Erosion 

Equation  1  is  typically  used  to  describe  surface  erosion.  In  this 
equation  there  arc  two  parameters  that  are  not  easily  defined:  the 
erosion  rate  coefficient  (M),  and  the  critical  shear  stress  of  erosion 
initiation  (x  ).  Both  have  been  predicted,  using  various  methods,  in 
the  laboratory.  A  field  experiment  was  performed  to  evaluate  both  the 
applicability  of  Equation  1  and  the  usefulness  of  the  laboratory  pre¬ 
dictions  of  these  parameters.  In  support  of  Equation  1,  we  find  that  a 
linear  equation  best  relates  the  erosion  rate  to  the  bed  shear  stress. 

A  comparison  was  made  of  the  measured  values  of  both  the  critical  shear 
stress  and  the  erosion  rate  coefficient  to  those  values  predicted  in  the 
laboratory.  We  found  that  flume  experiments  using  naturally  deposited 
sediment  beds  do  a  reasonable  job  in  predicting  these  parameters,  but 
that  experiments  using  remolded  sediments  do  not  approximate  the  condi¬ 
tions  encountered  in  the  field. 

Deposition 

Equation  2  is  typically  used  to  describe  the  deposition  of  cohesive 
sediments.  The  parameter  described  as  the  critical  shear  stress  of 
deposition  is  rather  ill-defined.  It  refers  to  a  bed  shear  stress  at 
which  all  of  the  suspended  sediment  load  will  eventually  fall  out  of 
suspension.  All  of  the  values  presented  in  Table  2  were  obtained  in 


flumes  that  maintained  a  constant  shear  stress.  Due  to  the  nature  of 
field  experiments  (i.e.,  no  control  over  the  shear  stresses)  it  is  not 
possible  to  measure  this  parameter  directly. 

Transportation 

Equation  5  is  typically  used  to  describe  the  t ransportat ion  of 
suspended  sediments.  Modeling  difficulties  are  encountered  when  one 
attempts  to  define  the  vertical  distribution  of  the  suspended  sediment 
concentration  using  Equation  6.  Two  parameters  in  this  equation  pose 
difficulties:  the  settling  velocity  and  the  reference  concentration. 

The  settling  velocity  is  typically  assumed  to  be  a  constant  for 
modeling  purposes.  This  is  not  in  keeping  with  the  fluctuations  seen  in 
the  apparent  settling  velocity  data  that  were  calculated  from  the  field 
data.  The  apparent  settling  velocity  changed  more  than  one  order  of 
magnitude  over  the  duration  of  a  tidal  cycle.  The  assumption  that  the 
settling  velocity  is  a  constant  is  certainly  an  oversimplification.  We 
have  found  no  viable  means  of  predicting  the  settling  velocity  in  the 
laboratory.  The  standard  means  of  estimating  the  "constant  settling 
velocity"  in  the  laboratory  (still-water  settling)  does  not  take  into 
account  internal  shearing,  which  is  the  dominant  process  in  determining 
the  field  settling  velocities  of  flocculated  materials. 

Although  the  issue  of  reference  concentration  was  not  directly 
addressed  in  this  experiment,  the  following  observation  can  be  made 
based  on  the  data  collected.  It  appears  that  there  is  some  background 
level  of  suspended  sediment  concentration  that  was  returned  to  after 
each  resuspension  event.  Since  most  resuspension  events  do  not  involve 
the  entire  water  column  (at  least  in  artificially-deepened  estuaries), 
there  must  be  some  level,  close  to  the  surface,  that  maintains  a  rela¬ 
tively  constant  level  of  suspended  sediment  concentration.  It  is  also 
likely  that  the  concentration  at  this  level  would  be  reasonably  uniform 
over  a  wide  expanse  of  the  estuary.  This  is  just  speculation,  but  it 
suggests  that  there  may  be  a  possibility  of  obtaining  the  reference  con¬ 
centration  with  a  single  set  of  measurements  specifying  the  background 
concentrat ion . 


CONCLUSIONS 

1.  A  linear  equation  best  relates  the  erosion  rate  to  the  bed  shear 
stress : 


2.  Flume  experiments  using  naturally  deposited  sediment  beds  do  a 
reasonable  job  in  predicting  the  values  of  both  the  critical  shear  stress 
and  the  erosion  rate  coefficient.  Laboratory  test  using  remolded  sedi¬ 
ments  do  a  poor  job  of  predicting  these  values. 

3.  The  settling  velocity  is  normally  assumed  to  be  a  constant,  but  this 
research  suggests  that  the  settling  velocity  can  vary  by  over  an  order 
of  magnitude  during  a  single  tide  cycle. 
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Table  1.  Soil  Shear  Strength  Zones  as  Measured 
in  the  Laboratory  (Parchure  and  Mehta, 
1985) 


Shear  Strength 

Range 

Zone 

(dynes/cm  ) 

T 

0.4  to  1.7 

SO 

I 

2.1  to  6.2 

sc 

X 

5.7  to  6.7 

sm 

Table  2.  Cited  Values  for  Critical  Shear 
Stress  of  Deposition 


Reference 

Tcd 

Krone  (  1962 ) 

Ariathurai  ( 1978) 

Mehta  (1986) 

(dynes/c^ 

0.6 

0.2 

1.5 

1.0 

0.8 
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Figure  4.  Tide  data  for  the  duration  of  the  experiment. 


Figure  7.  Stations  at  which  measurements  were  recorded  during 
profile. 
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gure  12.  Apparent  settling  velocity  versus  a  referenced  sediment 
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Appendix  A 
VELOCITY  DATA 


This  appendix  contains  a  graphical  presentation  of  the  time  series 
velocity  data.  The  data  were  presen;  d  as  two  horizontal  velocity 
vectors  at  the  indicated  heights  (16  and  62  cm)  above  the  z-datum. 

The  time  scale  runs  continuously  from  Figure  A-l  through  Figure 
A-5.  The  time  0000  hours  on  Figure  A-l  is  midnight  at  the  beginning  of 
the  first  day  of  the  measurements.  The  time  5500  hours  on  Figure  A-5 
(for  example)  is  55  hours  later  (i.e.,  7:00  a.m.  the  third  day). 

An  averaging-type  curve  smoothing  routine  was  used  in  which  each 
velocity  value  was  modified  by  the  two  adjacent  values.  The  lines  shown 
are  smooth  curves  drawn  through  the  given  data  points. 

Vertical  velocity  distributions  for  the  profile  data  can  be  obtained 
by  using  Equation  7  in  the  main  text  of  this  report  with  the  friction 
velocity  and  roughness  coefficient  values  given  in  Appendix  B. 
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Appendix  B 


FRICTION  VELOCITY 


This  appendix  contains  a  graphical  presentation  of  the  time  series 
and  profile  friction  velocity  data  (Figures  B-l  through  B-5)  as  well  as 
a  list  of  the  profile  data  (Table  B-l).  The  profile  data  list  (Table 

B-l)  contains  the  friction  velocity  and  the  roughness  coefficient,  both 

of  which  are  based  on  the  measured  vertical  distribution  of  the  horizontal 
velocity  vector. 

The  friction  velocity  time  series  data  were  computed,  using  Equation 

7  in  the  main  text  of  this  report,  from  the  velocity  time  series  data 

shown  in  Appendix  A.  The  time  scale  runs  continuously  from  Figure  A-l 
through  Figure  A-5.  The  time  0000  hours  on  Figure  A-l  is  midnight  at 
the  beginning  of  the  first  day  of  the  measurements.  The  time  5500  hours 
on  Figure  A-5  (for  example)  is  55  hours  later  (i.e.,  7:00  a.m.  for  the 
third  day).  The  line  shown  is  a  smooth  curve  drawn  through  the  given 
data  points. 
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Table  B-l.  PROFILE  DATA 


TABLE  OF  FRICTION  VELOCITIES  AND  ROUGHNESS  COEFFICIENTS 


PROFILE 

FRICTION  VELOCITY 

ROUGHNESS  COEFFICIENT 

( day -hr  ) 

( cm/sec  ) 

( cm/sec  ) 

266-1530 

1  .55 

0.47 

266-1600 

1.95 

0.49 

266-1630 

1  .24 

0 

266-1700 

1  .75 

0.002 

266-1730 

0 .67 

0 

266-1800 

1  .52 

0.11 

267-0700 

1  .  00 

0.19 

267-0730 

0 .64 

0 

267-0830 

1  .47 

6.83 

267-0900 

0.12 

0 

267-0930 

0.83 

0.03 

267-1000 

1 .37 

0.31 

267-1030 

0.74 

0 

267-1100 

0.91 

0.05 

267-1130 

1  .02 

0.18 

267-1200 

0.82 

0.02 

267-1230 

0.66 

0.003 

267-1300 

0.50 

0.03 

267-1330 

0.58 

0.24 

267-1400 

0.83 

0.40 

267-1430 

1 .48 

8.60 

267-1500 

0.64 

0.46 

267-1530 

0.24 

0 

267-1600 

1  .23 

0.50 

267-1630 

2.97 

0.22 

268-0730 

0.43 

0 

268-0800 

0.40 

0 

268-0900 

1  .30 

4.66 

268-0930 

0.94 

1 .02 

268-1000 

0.88 

0.02 

268-1030 

1.71 

0.42 

268-1100 

1  .10 

0.03 

268-1130 

0.68 

0.06 

268-1200 

0.29 

0 

268-1230 

1  .10 

0.21 

268-1300 

0.30 

0 

268-1330 

0.81 

3  .38 

B-2 


FRICTION  VELOCITY 


B-3 


1500  1600  1700  1800  1900  2000  2100  2200  2300  2400 

TIME  (HRS) 

Figure  B-l 


B-4 


Csi 


400  2500  2600  2700  2800  2900  3000  3100  3200  3300 

TIME  (HRS) 

Figure  B-2 


B-5 


3200  3300  3400  3500  3600  3700  3800  3900  4000  4100 

TIME  (HRS) 

Figure  B-3 


B-6 


TIME  (HRS) 

Figure  B-4 


B-7 


5000  5100  5200  5300  5400  5500  5600  5700  5800  5900 

TIME  (HRS) 


Appendix  C 

SUSPENDED  SEDIMENT  CONCENTRATIONS 


This  appendix  contains  a  graphical  presentation  of  the  suspended 
sediment  concentration  time  series  data  (Figures  C-l  through  C-5).  The 
data  presented  are  that  of  five  scalar  concentration  values  measured  at 
the  indicated  heights  (13,  21,  34,  58,  and  87  cm)  above  the  z-datum. 
Frequently,  the  sensors  at  the  13  and  21  cm  levels  were  outside  their 
operating  range;  therefore,  these  points  are  only  shown  during  the  times 
of  operation. 

The  time  scale  runs  continuously  from  Figure  A-l  through  Figure  A-5 
in  Appendix  A.  The  time  0000  hours  on  Figure  A-l  is  midnight  at  the 
beginning  of  the  first  day  of  the  measurements.  The  time  5500  hours  on 
Figure  A-5  (for  example)  is  55  hours  later  (i.e.,  7:00  a.m.  on  the  third 
day).  The  lines  shown  are  smooth  curves  drawn  through  the  given  data 
points . 

Vertical  suspended  sediment  concentration  distributions  for  the 
profile  data  can  be  obtained  by  using  Equation  6  in  the  main  text  of 
this  report,  with  the  reference  concentration  and  settling  velocity  data 
presented  in  Appendix  D  and  the  friction  velocity  data  presented  in 
Appendix  B. 
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Appendix  D 

APPARENT  SETTLING  VELOCITY 


This  appendix  contains  a  graphical  presentation  of  the  time  series 
and  profile  settling  velocity  data  (Figures  D-l  through  D-5)  as  well  as 
a  list  of  the  profile  data  (Table  D-l).  The  profile  data  list  contains 
’jo ili  the  apparent  settling  velocity  that  is  based  on  the  measured 
vertical  suspended  sediment  concentration  distribution,  and  the  refer¬ 
ence  concentration  used.  The  reference  concentration  is  the  measured 
value  at  82  cm  above  the  z-datum.  Data  are  not  presented  for  the  latter 
portion  of  day  268  due  to  the  failure  of  the  concentration  sensor  at  the 
level  z  =  34  cm.  This  failure,  combined  with  the  fact  that  sensors  at  z 
=  13  and  21  cm  were  operating  outside  of  their  range,  reduced  the 
vertical  extent  of  the  profile  sufficiently  to  render  it  useless. 

The  settling  velocity  time  series  was  computed,  using  Equation  6  in 
the  main  text  of  this  report,  from  the  suspended  sediment  concentration 
data  presented  in  Appendix  C.  The  time  scale  runs  continuously  from 
Figure  A-l  through  Figure  A-5  in  Appendix  A.  The  time  0000  hours  on 
Figure  A-l  is  midnight  at  the  beginning  of  the  first  day  of  the  mea¬ 
surements.  The  time  5500  hours  on  Figure  A-5  (for  example)  is  55  hours 
later  (i.e.,  7:00  a.m.  on  the  third  day).  The  line  shown  in  the  graphical 
presentation  is  a  smooth  curve  drawn  through  the  given  data  points. 
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Table  D-l.  PROFILE  DATA 


PROFILE  DATA 


TABLE  OF  SETTLING  VELOCITIES  AND  REFERENCE  CONCENTRATIONS 


PROFILE 

SETTLING  VELOCITY 

REFERENCE  CONCENTRATION 

( day -hr  ) 

( cm/sec ) 

( mg/ 1  ) 

266-1530 

.067 

68.4 

266-1600 

.  069 

81  .1 

266-1630 

.163 

116.3 

266-1700 

.039 

73.2 

266-1730 

.001 

87.1 

266-1800 

.026 

94 .9 

267-0700 

.016 

75 . 7 

267-0730 

.012 

73  .9 

267-0830 

.016 

71  .9 

267-0900 

.003 

69 .1 

267-0930 

.024 

64.7 

267-1000 

.034 

68.4 

267-1030 

.020 

69.6 

267-1100 

.031 

72.5 

267-1130 

.014 

77.5 

267-1200 

.058 

227.9 

267-1230 

.069 

145.9 

267-1300 

.048 

144.2 

267-1330 

.011 

97.1 

267-1400 

.009 

84.5 

267-1430 

.052 

92.9 

267-1500 

.056 

126.5 

267-1530 

.027 

106.3 

267-1600 

.275 

74.3 

267-1630 

.337 

114.2 

268-0730 

.011 

73.5 

268-0800 

.003 

69 .6 
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